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Abstract. In this paper, we introduce a simplified model for explain- 
ing the observations of the optical intraday variability (IDV) of the BL 
Lac Objects. We assume that the source of the IDV are the stochas- 
tic oscillations of an accretion disk around a supermassive black hole. 
The Stochastic Fluctuations on the vertical direction of the accretion 
disk are described by using a Langevin type equation with a damping 
term and a random, white noise type force. Furthermore, the prelim- 
inary numerical simulation results are presented, which are based on 
the numerical analysis of the Langevin stochastic differential equation. 
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1. Introduction 

Intraday variability is usually defined as the variation of the differential photometric 
magnitude greater than 3 times of sigma value. However, if the variation is greater 
than 5 times of sigma value is the best and reliable (Wagner and Witzel 1995 ; 
Fan, J. H. 2005 ; Gupta et al. 2008). There are many observations for the intraday 
variation in difference bands, and to explain the observational data many theoretical 
models have also been proposed for explaining these interesting phenomena. A 
cellular automaton model was introduced initially by Kawaguchi et al. (1998), and 
further developed in Tang (1999). 

It is the purpose of the present paper to propose an alternative model for the 
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explanation of the observed intra-day variability in BL Lac objects and for other 
similar transient events. The basic physical idea of our model is that the source of 
the intraday variability can be related to some stochastic oscillations of the disk, 
triggered by the interaction of the disk with the central supermassive black hole, 
as well as with a background cosmic environment, which perturbs the disk. To 
explain the observed light curve behavior, we develop a model for the stochastic 
oscillations of the disk, by taking into account the gravitational interaction with the 
central object, the viscous type damping forces generated in the disk, and a stochas- 
tic component which describes the interaction with the cosmic environment. The 
mathematical model is formulated in terms of a stochastic, Langevin type differ- 
ential equation, describing the stochastic oscillation of the accretion disk. This 
stochastic oscillation model can reproduce the aperiodic light curves associated 
with transient astronomical phenomena. Usually the IDV phenomena is explained 
by using the jet model. There are definitely interactions between the jet and the 
accretion disk. The stochastic fluctuations of the disk can influence the jet through 
some energy transfer, and they can represent a source for the jet perturbations. 

2. Vertical thin disk undamped oscillations 

We consider the oscillations of the disk as a whole body under the influence of 
gravity. The disk is considered thin in the sense discussed by Shakura and Sun- 
yaev (1973). We assume that at the center of the disk we have a compact object 

of mass M. We approximate the distribution of the surface density in the geo- 
metrically thin disk by the formulae S = Sq = constant, Rin < r < Radj^ and 



S = Sq ("R^j , Radj < f < Rout, respectively, (Shakura and Sunyaev 1973; 



Shakura and Sunyaev 1976; Titarchuk and Osherovich 2000), where Rin is the in- 
nermost radius of the disk, Radj is an adjustment radius in the disk, and i?out is the 
outer radius of the disk. The index 7 of the surface density can be either 3/5 or 



The vertical oscillations of the disk as a whole can be described by the equation 
of motion M^d^z/dt^ + Fo{z) = 0, or, equivalently. 
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3. Stochastic fluctuations of accretion disks 



In the following, we consider a model in which we assume that the accretion disk 
around a massive central object is immersed in a fluctuating fluid flow. For sim- 
plicity describe the accretion disk as a macroscopic sphere of radius a. In order 
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to describe the stochastic processes in oscillating disks, we replace the equation of 
motion of the disk, given by Eq. (1), with a stochastic Langevin type differential 
equation, which can be written as 

(Pz dz 2 Fit) 

where the term ^dz / dt, ^ =constant, takes into account viscous dissipation. The 
fluctuating force F{t) is independent of z, and it varies extremely rapidly as com- 
pared to z. Since F{t) is very irregular, we assume that its mean value is zero, 

im) = 0. 

By defining 

as the total energy of the oscillating disk, the luminosity of the disk, representing 
the energy lost by the disk due to viscous dissipation and to the presence of the 
random force, is given by 

dE ^ ^ rrfdzV dF 

As an example of application of this model to a concrete astrophysical sys- 
tem let's consider the case of a thin disk around a compact general relativistic ob- 
ject of mass M = IO^Mq. We assume that the adjustment radius of the disk is 
Radj = 75Rin, while the outer radius of the disk is located at Rout = 250i?j„. By 
assuming 7 = 3/5 we obtain first from Eq. (2) the frequency of oscillations of the 
unperturbed disk as loq = 1.89 x 10~^ s~^. The corresponding oscillation period 
of the disk in the absence of any dissipative or random processes is T = 18.739 
hours. For the mass of the disk we take Md = 1.58 x 10^^ g = 7.9Mq. For the 
dumping coefficient ^ we obtain = 6TrHr]/Md = 2.65xl0~^s~^, where for the 
viscosity coefficient we have assumed the value rj = 10^^ erg s/cm^. Therefore the 
equation of motion of the stochastically oscillating disk is given by 

z dz 
—r + 2.65 X 10"^— + 1.89 X IQ-^z = 6.329 x W-^^F(t). (6) 
at^ dt 

The results of the numerical integration of Eq. (3) are represented in Figs. 1 and 



4. Conclusions 



We have developed a theoretical model describing the stochastic oscillations of the 
disk. A comparison of the model predictions with the observational data of IDV 
from different sources will be done in a forthcoming paper. For the case of the 
ADAF disks, we can also apply the Langevin equation in the vertical direction. 
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Figure 1. The stochastic force F{t)/M4 cm/s^ (left figure) and z (right figure) as a function 
of time. 
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Figure 2. The stochastic velocity (left figure) and the luminosity of the disk (right figure) 
as a function of time. 

since this type of disks also has some vertical oscillation modes. Although the 
geometry of the ADAF disk is geometrically thick and optically thin, it can also 
provide a source for the disturbances of the jet. The effect of the jets in the frame- 
work of the present model and the effects of the relativistic corrections to the model 
will be also analyzed. 
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